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Abstract 

"We derive the distribution of energy and momentum transmitted from a primary fast parton and its medium-induced bremsstrahlung 
7-H gluons to a thermalized quark-gluon plasma. Our calculation takes into account the important and thus far neglected effects of 
CO quantum interference between the resulting color currents. We use our result to obtain the rate at which energy is absorbed by the 

medium as a function of time and find that the rate is modified by the quantum interference between the primary parton and sec- 
>^ondary gluons. This Landau-Pomeranchuk-Migdal type interference persists for time scales relevant to heavy ion phenomenology. 

We further couple the newly derived source of energy and momentum deposition to linearized hydrodynamics to obtain the bulk 
^> jnedium response to realistic parton propagation and splitting in the quark-gluon plasma. We find that because of the characteristic 

large angle in-medium gluon emission and the multiple sources of energy deposition in a parton shower, formation of well defined 

Mach cones by energetic jets in heavy ion reactions is not likely. 
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1. Introduction 

The suppression in the production rates of energetic leading 
particles and jets in relativistic heavy ion reactions relative to a 
naive superposition of nucleon-nucleon collisions is one of the 
most striking results from the heavy ion program at the Rela- 
tivistic Heavy Ion Collider (RHIC) |Q]] and now at the Large 
Hadron Collider (LHC) JUS]. These phenomena, often collec- 
tively referred to as "jet quenching" |4[, have been studied ex- 
tensively both experimentally and theoretically and are thought 
to provide valuable information about the quark-gluon plasma 
(QGP) created in these events 

At a fundamental level, the physics of jet quenching is largely 
reflective of the interaction of energetic, or fast, partons with 
the QGP medium which they traverse. Fast partons lose energy 
primarily through medium-induced radiation fl. Thus, the 
original fast parton evolves into an in-medium parton shower. 
The shape and energy distribution associated with this shower 
provide more information about the underlying QCD dynamics 
than the suppression of leading particles alone and form the ba- 
sis for using full jet observables as a new and powerful probe 
of the QGP 0]. 

Parton showers can also be substantially modified through 
collisional energy losses to the underlying medium. Although 
the collisional energy loss associated with a single parton is 
thought to be relatively small, the cumulative effect associated 
with a full parton shower can become quite large as the radiated 
gluons themselves become sources of collisional energy loss. 
Thus, the question of how much energy a parton shower trans- 
mits to the medium is of fundamental interest in the description 
of full jet observables. 
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A closey related question is how the medium responds to a 
parton shower. As a parton shower loses energy and momen- 
tum, the underlying medium is affected in a way that depends 
on its properties as well as the space-time distribution of the 
energy and momentum deposition. This problem has gained 
attention due to experimental measurements of azimuthal par- 
ticle correlations associated with high pj triggers in heavy ion 
collisions that display a double-peaked or conical structure Jit]. 
These measurements may reflect the interaction of fast par- 
tons with the medium and among the proposed explanations 
for the structure are Mach cone Shockwaves excited in the bulk 
medium by fast partons (2D - Other explanations that do not re- 
flect the interaction of fast partons with the medium, such as 
fluctuating initial conditions and triangular flow have also been 
proposed [10]. Whether or not the conical structure associated 
with high pr triggers at RHIC is associated with the medium 
response to a parton shower, this topic remains important, par- 
ticularly in the light of the the new exciting possibilities for jet 
physics with the LHC heavy ion program. 

Both the question of how much energy a parton shower trans- 
mits to the medium and how the medium responds to a parton 
shower can be addressed in the same framework by calculating 
the source term associated with the parton shower. The source 
term, here denoted J v (x), is the space-time distribution of en- 
ergy and momentum flowing between a parton shower and the 
underlying medium. It couples to the energy-momentum tensor 
(EMT) r /JV as <9 A ,r /JV = J v . The source term carries information 
about the rate of energy transfer to the medium and acts as a 
source (hence the name "source term") for the evolution of the 
underlying medium in the presence of the parton shower. 

Previous calculations of the source term have focused on a 
single primary parton 111 ill . Applications to Shockwave for- 
mation that have attempted to include medium-induced gluon 
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bremsstrahlung have neglected the spatial extent of the source 
and, just as importantly, the color quantum interference be- 
tween the primary parton and associated radiation II12I . In what 
follows we will perform the first calculation of the source term 
associated with a parton shower. Our calculation will include 
the full spatial extent and color quantum interference effects, 
and will be presented in closed analytic form. To perform this 
calculation, we will start with the result of a recent paper by 
Neufeld 1.13(1 in which the source term induced by a single fast 
parton in a medium of thermal quarks was derived in integral 
form using thermal field theory. The source term was obtained 
by taking the thermal average of the divergence of the quark 
energy-momentum tensor, with the fast parton coupled to the 
medium by adding an interaction term to the Lagrangian. 

Our paper is organized as follows: we present the theoreti- 
cal formalism for the evaluation of the source term associated 
with a parton shower in section [2] In section [3] we present two 
applications of the source term we derive: the rate of energy 
transfer between the parton shower and the medium and the re- 
sponse of the medium to the parton shower. We show that a 
significant amount of energy can be transferred from the par- 
ton shower to the medium depending sensitively on the average 
angle of gluon radiation. We further show that for realistic av- 
erage angle of gluon emission the medium response no longer 
appears as a well defined Mach cone but rather a superposition 
of several distinct perturbations in the medium. Our summary 
and conclusions are given in section |4] 



2. Theoretical formalism 

Our approach in this section is to extend the result pre- 
sented in H 1311 to the case of an asymptotically propagating 
primary fast parton which undergoes medium-induced gluon 
bremsstrahlung. In this manuscript we use the following 4- 
vector notation: x = {t,r),u = (l,u). 

Our starting points are the Feynman diagrams of Figure [T] 
The dark blobs represent the divergence of the medium energy- 
momentum tensor, for which we will implement the result de- 
rived in 1U3I1 . In the left diagram of Figure Q] the gluon lines 
labeled with p4 and p$ connect the medium to the primary fast 
parton (a quark in the diagram), in the middle diagram they con- 
nect to the radiated gluon, and in the right diagram they connect 
to both. This diagram represents the non-trivial quantum inter- 
ference between the two partons. A fourth diagram with p\ and 
ps exchanged in the interference term is not shown in Figure [1] 
but is included in the calculation presented here. 

In Figure [T] the SU(3) color matrices are denoted by upper 
and lower case T and t for the adjoint and fundamental repre- 
sentations, respectively. Using standard color algebra, we find 
that the color factors for the first two diagrams in Figure Q~|sim- 
ply yield the quadratic Casimir of the particular representation: 
Cf = 4/3 for a quark (first diagram) and Ca — 3 for a gluon 
(middle diagram). The color factors for the third and fourth (not 
shown) diagrams yield —Ca/2 = -3/2 each. We point out that 
the color factors of these last two contributions remain the same 
if the primary parton is a gluon. Note that we have averaged 



over the initial colors and factored out the Casimir associated 
with the bremsstrahlung vertex. 

If for the moment we write the divergence of the medium 
EMT as an arbitrary function J v (x, 111,112) (with u\ and mt being 
the velocity 4-vectors of the primary parton and radiated gluon, 
respectively), then the combined result for the four contributing 
diagrams is: 

dfj^ - C p J v a (x, mi, mi) + C A J v (x, 1*2, ui) 

-^-[J y (x,M U u 2 ) + J v (x,U2,M l )] , (1) 

where C p is the quadratic Casimir of the primary parton and 
the subscript a indicates that we consider the primary parton to 
propagate asymptotically. In the limit of perfectly collinear ra- 
diation one has u\ = 112 and the interference terms exactly can- 
cel the contribution from the radiated gluon. This result reflects 
the fact that the medium will never resolve the bremsstrahlung 
gluon from the parent parton and the color of the system is the 
same before and after the radiation. However, if radiation is not 
perfectly collinear (and u\ + u%), the system evolves in time, 
the two new color currents separate, and the interference terms 
gradually vanish. At t — » +00 the system begins to look like 
two separate color charges. 

In this work we consider the primary parton to be propa- 
gating asymptotically, that is we ignore any finite time effects 
associated with the initial large Q 2 scattering in a heavy ion 
collision. This problem has been examined in a different con- 
text in 111 411 . Although these effects are important to understand 
and may very well be significant, our focus is instead on the 
time dependence associated with the gluon radiation and the 
interference effects that arise with it. When both gluon lines 
attach to the primary quark, the current can be represented by 
the classical asymptotic form: f 1 — > gQ a p i/ y 5(r - vl\ t) with 
Qp Q", - C p . Similarly, when both gluon lines attach to the 
secondary gluon the current can be represented by the form 
jf _> gQa M ^ gfj. _ U2 t y However, when the gluon lines 
are connected as in the third diagram of Figure [1] there is no 
simple way to express the current and one must combine the 
product. The result is reflected in Eq. ([T}: j M (x\) j v *(x2) — * 
-@(fi)®(?2)g 2 %w^(ri - ui ti)u\6(r 2 - \i2tz), where we only 
consider interference effects after the point of emission. 

With these considerations in place, we write the integral ex- 
pression for the source term derived in 111 311 for an asymptot- 
ically propagating fast parton in a medium of thermal quarks, 
corresponding to the first diagram of Figure Q~]and written in the 
notation of Eq. ([TJ: 



/;(*,«!,«!) = -4>W f dp3 *% dp5 e-*"+»>n F (p 

J (27T) 

x 6(pl)G R (p 4 )G R (p3 + p4)G R (p 5 )6(p 4 ■ ui)S(p 5 ■ u\) 



) 



[(2(/»3 • «l) 2 - u\p3 ■ P4)p V 5 
-U\(p3 ■ Wl)(2/? 3 ■ p 5 + p 4 ■ p 5 )j , 



(2) 



where Gr(p) = (p 2 + iep Q Y x is the retarded Green's function, 
rif(p) = (e' p ^ T + 1) is the Fermi distribution function, T is 
the temperature, and Nf is the number of active flavors. 
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Figure 1: Feynman diagrams for the source term induced by a primary parton that undergoes medium-induced bremsstrahlung. The color structure of the different 
contributions to the source term is explicitly shown. The third diagram represents the non-trivial interference between the parent parton and the radiated gluon (the 
second way of attaching p4 and p$ to the source is not shown explicitly). 



For our calculation we also need the generalization of Eq. (0 
to the case of two separate velocities, that is J v (x,u\,u 2 ), cre- 
ated at an initial time t — 0. The primary effect of this modifi- 
cation enters through the 5 function structure of Eq. (0: 



6(p4 ■ u x )5(p 5 - Mi) 



1 



(2n) 2 (p 4 ■ u 2 + ie)(p5 ■ u\ + ie) 



(3) 



The above replacement can easily be understood from the inte- 

poo 

gral representation of the 6 function: 2n5{p ■ u) = J_ ate p ". 
For the case of a charge created at time t — the integration be- 
gins at rather than -co. One can see that after contour integra- 
tion in the p° A , p° 5 variables the primary result of the replacement 
made in Eq. (0 is that p4 and p$ no longer enter the exponential 
symmetrically when u\ + u 2 . The result is that the contribution 
from J v (x, u\,u 2 ) will die out on an timescale related to the in- 
verse of the difference between u\ and u%. 

The specific generalization of the polynomial structure in the 
numerator (the third and fourth lines) of Eq. (0 to the case of 
two separate velocities is not as essential as the replacement 
shown in Eq. (or the pole structure). In particular, correc- 
tions to the form of Eq. are at most proportional to Ui - U2, 
which should be suppressed when the radiation angle is smaller 
than one. With this in mind, we make the obvious assignment 
that the first term remain symmetric in u\,u 2 . For the second 
term we use the form: 



u\(p3 ■ ui)6(p 4 ■ ui)6(p 5 ■ «i) 



-u[(p 3 ■ u 2 ) 



(2n) 2 (p4 ■ m + ie)(p 5 ■ u\ + ie) 

(4) 



based on the relation of pi and p4 in the Green's function struc- 
ture. In the worst case, the correction in this second term is 
proportional to Ui - 112. The combined result can be written as: 



T\x,u\,u 2 ) = 4iNFg 



d 4 p 3 d 4 p 4 d 4 p 5 _ ix . (m 



Mf3) 



8(pI)Gr(P4)Gr(P3 + P4)G R (p 5 ) 



(P4 ■ u 2 + ie)(p 5 ■ u\+ ie) 

X [(2(p 3 • Ml)(p 3 • U 2 ) -Ui-U 2 p3- Pa)P V 5 

-u\(p3 ■ u 2 )(2p 3 ■ p 5 + P4 ■ Ps)\ ■ (5) 

The full source term, which is the sum of the terms in Eq. dl}, 
can now be obtained by evaluating the expressions given in 



Eqs. (f2]i and (0. As we will show in what follows, one ac- 
tually only needs to evaluate Eq. (0, for which Eq. (0 can be 
obtained as a special case. 

We now proceed to evaluate Eq. (0. In order to simplify the 
calculation we use bare propagators and take the hard thermal 
loop (HTL) approximation. These simplifications will allow us 
to obtain an analytic result for the source term. We also will 
consistently ignore radiative effects arising from the formation 
of a charge at an initial time. These effects are associated with 
causality and are in addition to the medium induced radiation 
we are considering in this paper. We postpone their considera- 
tion to a future work. We emphasize that all vacuum contribu- 
tions have already been subtracted from Eq. (0. Furthermore, 
although 01 was derived for a medium of quarks (as can be seen 
from the number of active flavors Nf and the Fermi distribution 
rif(p)) the extension to include medium gluons is trivial in the 
HTL limit as will be shown below. 

We obtain the HTL approximation by taking the limit that 
the fields generated by the external current are soft compared 
to the temperature, or formally \p 4 \ ~ gT <K T in Eq. 0. The 
coupling constant g is formally considered to be much less than 
one. Explicitly, we expand Gr(/73+/?4) in the limit of |p4| <s \pj,\ 
since \p$\, which appears in the Fermi distribution function, is 
cut off in the integration at a value on the order of T to low- 
est contributing order. We find that in the HTL approximation 
Eq. 01 reduces to: 

G R (p4)G R (p 5 ) 



(P4 ■ u 2 + ie)(p 5 ■ ui + ie) 
p\ ((v ■ «i)(v ■ u 2 ) p v 5 - u\ (v ■ u 2 )(v ■ p 5 )j 



(P4-V + ie) 2 
ui ■ u 2 v ■ p4 p v 5 + u\ (v ■ u 2 )(p 4 ■ ps) 



P4 ■ v + ie 



(6) 



where m 2 D = g 2 T 2 Nf/6 for a thermal medium comprised 
quarks only. In Eq. (0 dQ. is the integration measure over 
the solid angle defined by the unit vector v, and v* 1 = (l,v). 
As mentioned previously, the extension to include medium 
gluons is straightforward and is done with the modification 
ml^g 2 T 2 (l +N F /6). 



In order to evaluate Eq. © it is convenient to note that the p4 
and integrations separate. The p$ contribution has the same 
general form in each term which we evaluate as: 



(In) 4 (p 5 ■ m + ie) 



Gr(ps) 



-0(f 2 



2 n (ti(«i ■ r " «T 0. r - ui y?(* - «i ■ r)) 
r ) An(r 2 -t 2 + y\(t - ui ■ r) 2 ) 3 / 2 ' (?) 



where j\ = 1/ ^1 - Wj, ? > 0, and we ignore radiative effects 
on the light cone associated with the formation of the field. The 
result of Eq. (O can be derived in a straightforward manner by 
expressing components of ps as a derivative on x and using the 
relationship that J d 4 pe~ lx ' p G R (p) = (2^) 3 5(f 2 - r 2 ) along with 

(p 5 • Mi + ze) _1 = -i C° dwe' wps ' Ul . The step function in Eq. © 
is a result of causality and ensures that the final result for the 
source term will only contribute in the time-like region. 

The p4 integration in Eq. (0 cannot be separated from the v 
integration and thus is more involved. For the first term one is 
aided by the cancellation occurring from ph Gr(pa) = 1. Again 
using (f>4 • u 2 + ie) - —i : L dwe' wpi '" 2 and writing (p4 ■ v + 
ie)" 2 - -d I dp%p4 - v + ie)~ x followed by an integration by parts 



allows: 



1 



g(y) 



(P4 ■ u 2 + ie) (p4-v + ie) 2 

i,Hix - U2K - VT \t - w) g(v) , (8) 



r d 4 p 4 r do. . 

J (In) 4 J ~An~ £ 

— i I dwdr - — 

Jo J (2^) 4 J 4/r 

where g(v) = ((v • u\){v ■ u 2 ) p v 5 - u\ (v ■ u 2 )(v ■ ps)j- At this 
point the p4 integration is performed to yield a 6 function, 
which enables the remaining integrals to be performed analyti- 
cally. The result is obtained as 

J (2n) 4 J 4tt e ( P 4 ■ u 2 + ie) (p 4 • v + ie) 2 
iy 2 g(v = v')0(t 2 -r 2 ) 



An 



yjr 2 -t 2 + y\(t - u 2 • r) 2 



(9) 



where 
i 

v = 



r - u 2 f 



{y 2 (xi 2 -r-u 2 t) + y 2 ^r 2 - t 2 + y\(t - u 2 • r) 2 j 



+ u 2 



(10) 



Lastly, we must evaluate: 
j4 - r dQ.e- ixPi G R (p 4 ) 



C <f PA r 

J (2tt)8J 



47r (/?4 • u 2 + /e) 



mi ■ u 2 p 5 + 



u v Av ■ u 2 )(p4 ■ Ps) 



P4 ■ v + ie 



(11) 



The first term in Eq. ([TP can be obtained in the same manner 
as Eq. (0 with the result: 



d 4 P 4 e- ixp *G R ( P 4) 
(2n) 4 (p4 ■ u 2 + ie) 



iy 2 0(f 2 - r 2 ) 



An Jr 2 -t 2 + y\(t - u 2 ■ r) 2 



(12) 



The second term is more involved and for the sake of brevity 
we briefly outline the important steps in obtaining the result. 
We first write (p4 ■ v + ie)~ l = -iC dTe' TPi ' v and then use 
contour methods to perform the p9 integration (which places a 
restriction that r < t). One can next perform the r integration 
to leave: 



f d 4 P 4 f 
J (2n) 4 J 

/dpd: 



dQ e- ixp >G R { P 4) 
An (p4 ■ u 2 + ie) 

dQ.i dQ. 



(v • u 2 )p a A 



v{u 2 ■ v)- 



1 



) 3 An 



„itp\y(u 2 -\) 



Im 



iv a e -ip(x-V!) 



P4 -v + i e 

(1 _ gitpCvi-v)) 



(U 2 -Vi) (vi-v) 



(u 2 ■ v) 



V] • (u 2 - v) 



-e 



ip\[-(x-u 2 t) 



(U 2 ■ Vi,Vi)" 



(U 2 -Vl) 2 -1 



(13) 



The second line of Eq. (fTTt can be simplified by redefining 
the v angular variable in terms of vi ■ v: that is, let Vi ■ v = cos 9, 
where 6 is associated with dQ. One can further use the rela- 
tions J dQ.e ,kr = Ansinkr/kr and J dD.ke' kr = Anir(sinkr - 
kr cos kr)/kr 3 to reduce the expression to two integrals which 
can be performed in closed form using standard integral rela- 
tions. The result is: 



/dp d: 



dCli dQ. i(l,\i)e- ip ^ Vl) (1 - e itp(vvv) ) 

Im : ; ; ; (u 2 • v) 



) 3 An' 



(« 2 ' Vl) 



— ; 
Anr 



(vi • v) 

-1 rri 



V-/tanh -1 [£]^ r fr - f tanh _1 [f] 



(14) 



where for simplicity we have restricted our attention to the re- 
gion of r < t knowing from Eq. (0 that this is the only region 
that contributes in the final expression. 

Finally we consider the third line of Eq. ([T3V which can 
again be simplified by redefining the v angular variable in 
terms of Vi ■ v. One can further use the integral relations 
Jo d4)\\, cos <fr\ exp ikb cos[4> — a] = 2n[Jo(kb), i xos aJ\ (kb)] 
and 11511: 



f 



dke iky Ji(kb) = f 1 1 + 



dke iky J (kb) 



1 



1 



yjb 2 -(y + ie) 2 

iy 



y/b 2 -(y + ie) 2 ) 



(15) 
(16) 



to again reduce the expression to two integrals. In the above 
expressions 7, are Bessel functions of the first kind, b > and 
e is a positive infinitesimal quantity. Performing the final two 
integrations, we find 



/dpd 



dQ 1 dQ / 1 
(u 2 ■ v) 



,!«pVr(U2-V) 



-e 



//>Vi-(X-U 2 



= l- 



2 2 
2 ' 



) 3 An " Vi ■ (u 2 - v) 

(u 2 • r) 2 , u 2 (r 2 - u 2 ■ r t) - r(u 2 ■ r ■ 



(U 2 ■ Vi.Vl) 



(u 2 -v,) 2 -l 
u\t)) 



-l- 



An 

u 2 - u 2 



(r - u 2 f) 2 r 2 - t 2 + yUt - u 2 ■ r) 2 



tanh" 1 - - / 



Antu 2 



u 2J 



+ i 



u 2 ■ r - uit 



u 2 f 



An{r-u 2 t) 2 47r(r-u 2 f) 2 , 



(17) 



We now have all of the necessary pieces to yield a final result 
for Eq. ((S). In order to achieve a compact form, is it helpful to 



define the notation w" — x a - y 2 {x ■ u^u", in which case one 
finds: 



Jo 



int. term: ti = 0.5 h = 2.5 [I'm] 



J v (x,u u u 2 )htl 



8/r 



j w\ i^5u ■ w 2 + u\ ■ u 2 j2 y]~ w ij 

( w +72 ^2j 2 ( _ w 2 )3 /2 



u\y\ y\{t-m-r) 



1 • "27271 w[ 



(-W 2 ) 3 / 2 



(_ w 2)3/2 



xdu 



(r - u 2 f) + 



(r - u 2 f) 2 

72 (u 2 (r 2 - u 2 ■ rf) + r(M 2 r - u 2 ■ r)) 



u 7i 



/ m 2 / _ u 2 ■ r + £u ■ u 2 y 2 (f - Ui • r) 



r — f tanlf 



tanh 1 m 2 - m 2 



'7?(ui ■ r - u\ t) (r • Ui)t?(/ - ui • r) 

1 1 + -z 



(18) 

where we have defined 6u = u\ — m 2 . Equation (fT8l is our 
primary result and we now discuss its general features and also 
how to specialize it to the asymptotic case J y a {x, u\,u\). Eq. ( fT8l 
is restricted by causality, which is contained in the step function 
0(f 2 - r 2 ). Furthermore, the first four lines of Eq. (fi"8) con- 
tribute in the asymptotic limit, whereas the last two lines are 
transient contributions, that is, they vanish as t — > oo even for 
Mi = m 2 . This can be seen by taking the case of u\ - m 2 and 
allowing u ■ r - u 2 t = const., while taking t — > oo. Thus one can 
obtain the asymptotic case J v a {x,u\,ui) by removing the step 
function and the final two lines of Eq. (fl8b . and of course set- 
ting Mi = m 2 . In this asymptotic limit our result should reduce 
to the result obtained previously by Neufeld in Bl ill . We have 
verified numerically that this is indeed true. 

Recall from Eq. ([TJ that the total source term is given by the 
combination C p J v a (x, u\,u\) + 3J v (x, m 2 , m 2 ) - j[J v (x, u\, ui) + 
J v (x, m 2 , Mi)], where the last two terms result from the interfer- 
ence graphs of Figure Q] and we refer to this contribution as 
"interference term" in what follows. To help illustrate the re- 
sult of Eq. dTSl we plot the interference term contribution for 
two different times in Figure [2] above. For all numerical re- 
sults we use the following parameters: g — 2, T — 0.35 GeV 
and mo = gT, which are based on the average values obtained 
for LHC collision energies in a Bjorken expanding plasma and 
we take Nf = for a gluon-dominated medium. Addition- 
ally we impose medium-induced Debye screening through an 
exponential decay factor, e~ mD(pi+P2 ^ 2 , where p\ and p 2 are dis- 
tances transverse to the axes of propagation defined by mi and 
m 2 , respectively. The velocity of the primary parton is chosen 
along the z axis with a magnitude determined by the energy and 
mass of the parton: E = ym with m — gT/ V6 for a quark and 
m — gT I y2 for a gluon. The magnitude of the velocity of the 




l [fmj 



Figure 2: A 3D representation of the extent and magnitude of the appropriately 
scaled interference term Jo/m D at two different times ; = 0.5 fm and t = 2.5 fin 
for a primary parton that propagates in the z direction. Note that this term 
disappears in the long time limit. 



radiated gluon is determined in the same way as the primary 
parton, and its direction determined by the angle of emission 
relative to the direction of the primary parton's propagation. For 
simplicity we will take the emission to occur in the x - z plane 
unless otherwise stated. Thus, all quantities are determined by 
specifying the energies of the primary parton and the radiated 
gluon, as well as the angle of emission. 

We plot the interference term contribution for two different 
times in Figure [2] above. For this plot we choose a primary 
quark with energy 5 Gev, a radiated gluon with energy 1 GeV 
and angle of emission 9 = 0.6 (angles will always be given 
in radians). Although these energies are phenomenologically 
small, they are convenient for the purposes of plotting. The plot 
shows the zero component of the interference term contribution 
scaled to a dimensionless ratio for times t = 0.5 fm and t = 2.5 
fm (in the same figure), where time is measured from the mo- 
ment of emission. The plot shows that at early times when the 
two partons are close together their interference contribution is 
important, but at later times as they become separated the in- 
terference contribution diminishes and eventually vanishes. In 
the next section we consider two applications of the source term 
derived in Eq. ( fT8l . 

3. Applications: differential energy transfer to the medium 
and linearized hydrodynamics 

The main theoretical result of this paper is Eq. ([T8V which 
combined with Eq. ((TJ gives the source term for a primary par- 
ton and a single radiated gluon as a function of time, including 
quantum color interference effects. In this section, we consider 
two applications of our result: the evaluation of the differential 
energy lost by a parton shower to the medium and the hydrody- 
namic response of the medium to a parton shower. 

The differential energy transfer to the medium is obtained 
from the source term through a spatial integration of the zero 
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Figure 3: Left panel: The rate of energy transfer from a q — > q + g system to the QGP medium for different gluon emission angles 8 = 0. 1, 0.3, 0.6. The naive 
no-interference result is presented by a solid line. Insert shows the total energy absorbed by the medium as a function of time. Right panel: The same rate and net 
energy transfer are presented for realistic multiple gluon emission of average angle 8 = 0.7 for both q — > q + 4 X g and g — » g + 8 X g systems. 



component: 



dE fro o 

— = I dr \C p Jq(x, mi, ui) + CaJ (x, ui, uq) 

- — (J°(x, M], U-l) + J°(X, M2, Ml)) 



(19) 



This transfer is of a collisional nature and the terms involv- 
ing only one velocity, C p J®{x, U\, u\) and CaJ°(x, 112,112), are 
fraught with an ultraviolet divergence typical of collisional en- 
ergy loss 111 611 . We regulate these terms with a standard short 
distance cutoff 0(p, - 1/(2 yJEfT)) related to the finite mo- 
mentum transfer, where p, is the distance transverse to the 
axis of propagation defined by m, and £, is the corresponding 
energy. On the other hand, the terms involving two veloc- 
ities, J°(x, Mi, M2) and 7°(x, Mi,M2), have no ultraviolet diver- 
gence and require no cutoff. However, motivated by the finite 
kinematics and for consistency, we apply a short distance cut- 
off to these terms reflective of two distinct gluon exchanges: 
0(pi - 1 /( V£7T))0(p2 - l/( yfE^T)). As these terms have no 
ultraviolet divergence, the sensitivity to the cutoff for the inter- 
ference contribution is small. 

We now consider two scenarios for the differential energy 
transfer to the medium which are both presented in Figure [3] 
The first is shown in the left panel and is for a primary quark of 
energy 75 GeV which emits a gluon of energy 4.5 GeV at time 
t — 1 .25 fm. At this point the primary quark's energy is reduced 
to 70.5 GeV. The differential energy loss to the QGP is shown 
for three different angles of emission, as well as without any 
interference effects. It should be noted that even in the absence 
of interference it takes on the order of y/mo for the radiative 
gluon to build up its asymptotic energy transfer rate (see the 
solid curve). The strength of interference effects depends on the 
transverse separation of the two color currents. When these cur- 
rents cannot be resolved by the medium the system behaves like 
the parent parton. This is illustrated by the smallest emission 
angle we consider, 9 = 0.1. The red dashed curve shows the 
fact that interference suppresses the energy loss to the medium 
for times extending up to about 3-4 fm after emission. For the 



larger angle of 9 — 0.3 there is a slight suppression initially fol- 
lowed by a slight enhancement persisting until about 2 fm after 
emission. Finally, we also consider 9 = 0.6 as we are moti- 
vated by the finding that the medium-induced bremsstrahlung 
is dominated by large angles 111 711 . For this largest angle the 
energy loss actually shows an enhancement relative to no in- 
terference which persists for about 2 fm after emission. The 
insert shows the total energy transferred to the medium for the 
different angles as a function of time. 

Next, we consider a more realistic representation of the 
medium-induced parton shower inspired by the Gyulassy- 
Levai-Vitev (GLV) energy loss formalism [18]. Such parton 
showering is relevant to phenomenological applications that 
aim to elucidate the medium response to jet propagation in 
hot and dense QCD matter II 1911 . Specifically, we make use 
of the bremsstrahlung spectra for quark and gluon jets aver- 
aged over the collision geometry in central Pb+Pb reactions 
at the LHC that have been previously employed to discuss jet 
and particle production in heavy ion collisions at the highest 
y/sNN lE3l- Medium-induced gluon emission is probabilistic 
and depends on the parton species. The most likely numbers 
are 4 (to g as 4.5 GeV) and 8 (io g * 4 GeV) for 75 GeV quark 
and gluon parent partons, respectively. We also note that the 
average emission angle is quite large, 6 as 0.7. The gluon emis- 
sion points are uniformly distributed along the parton propaga- 
tion path in a medium of length 5 fm. In the GLV approach 
the gluon distribution is azimuthally symmetric around the di- 
rection of jet propagation. For simplicity, we distribute them 
here in the x-z plane. Further, quantum interference effects are 
only taken between the primary parton and each radiated gluon 
- not between the multiple radiated gluons. The result is shown 
in the right panel of Figure [3] where again we find an enhance- 
ment of energy transfer to the medium. The insert shows the 
total energy transferred to the medium where one sees that for 
a gluon-induced shower its value is about 32 GeV for our pa- 
rameters. 

We now evaluate the hydrodynamic response of the medium 
to a parton shower. We consider the realistic scenario discussed 
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£e|x_| n n 

— ■ , - =0.08, 0=0.1, t = 5[fm] — ■ ' , - =0.08, 0=0.7, t = 5[fm] 
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Figure 4: In the left panel we show energy density disturbance generated by a parton shower originating from a primary quark with average gluon emission angle 
of 8 = 0.1, whereas in the right panel the average angle is 6 = 0.7. One can see that for the narrow emission angle the parton shower generates a well defined Mach 
cone similar to what one would expect from a single primary parton. However, the interference effects from such small angle emissions suppress the energy transfer 
to the medium. For the larger emission angle there is no longer a well defined Mach cone but rather a superposition of several distinct perturbations in the medium. 
The interference effects for this angle actually lead to an enhancement in the overall energy gained by the medium, even relative to no interference. 



directly above, only we now include both 6 - 0.7 and 6 = 0.1. 
By comparing these two angles we can illustrate clearly the ef- 
fect of the interference on the medium response. The imple- 
mentation of the full source term is rather numerically expen- 
sive therefore in this exploratory study we implement the sim- 
plest possible description of the shower by treating each par- 
ton as a 5 function source with a weight given by the energy 
loss rate to the medium: J v = [dEj/dt(t)] u^d{Tj - u,f), where 
the subscript ; refers to a summation over each parton in the 
shower. This simplified form of the source correctly accounts 
for energy and momentum conservation and allows us to in- 
corporate the energy loss rates evaluated above. We make the 
stipulation that the contribution coming from the quantum in- 
terference is associated with the radiated gluon rather than the 
primary quark. 

We focus on the energy density perturbation, 5e, generated in 
the medium by the parton shower. Within linearized hydrody- 
namics this quantity can be written in momentum space as [21] 



ikJ L (k) + J°(k)(ia) - T s k 2 ) 
a> 2 - c?k 2 + iT t wk 2 



(20) 



In the above equation, c s denotes the speed of sound and 
r, = jpf is the sound attenuation length with r//s being the 
shear viscosity to entropy density ratio in the medium. Also, the 
source vector has been divided into transverse and longitudinal 
parts: J = kJi + Jt- In addition to the parameters listed above 
we choose the proposed minimum bound for the shear viscos- 
ity to entropy density ratio rj/s = 0.08 12211 . Although shear 
viscosity estimates obtained from phenomenological studies of 
hydrodynamic simulations of the bulk matter produced in heavy 
ion experiments l23ll suggest a value perhaps 2-3 times larger 
than rjl s - 0.08, we have chosen the minimum bound because 



it is particularly well suited for plotting the important features 
of the medium response. 

In Figure [4] we present the result for 6e for the scenario 
and parameters discussed above. We consider a primary quark 
of energy 75 GeV which emits 4 gluons of energy 4.5 GeV. 
The gluon emissions are uniformly distributed over the primary 
quark's path length of 5 fm, at which point we plot the medium 
response. For the purpose of illustration we have chosen the 
gluon emissions to alternate between the positive and negative 
x directions, although in practice they are randomly distributed 
in the x - y plane. Additionally, we have scaled the energy 
density distributions associated with each parton by a factor of 
\x-\ = |r - ro - u(f - fo)| to compensate for the broadening of the 
Mach cones. Here, ro and to are the initial position and time of 
the particular parton. 

In the left and right panels of Figure|4]we show the result for 
6 = 0.1 and 6 = 0.7, respectively. One can see that for the nar- 
row emission angle the parton shower generates a well defined 
Mach cone similar to what one would expect from a single pri- 
mary parton ]24j]. However, the interference effects from such 
small angle emissions suppress the energy transfer such that the 
total energy gained by the medium is only moderately enhanced 
compared to the single primary parton (about 1.7 times larger 
for our parameters). Contrast this to previous studies which 
ignored the interference suppression and found significant en- 
hancement even for exactly collinear radiation [12]. For the 
larger emission angle there is no longer a well defined Mach 
cone but rather a superposition of several distinct perturbations 
in the medium. The interference effects for this angle actu- 
ally lead to an enhancement in the overall energy gained by 
the medium, even relative to no interference. Compared to the 
single primary parton the energy gain is about 2.7 times larger. 
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However, the price for this gain in energy is that there is no 
longer a well defined Mach cone structure. 

4. Conclusions 

In this Letter we evaluated the source term J v of energy and 
momentum deposited by a parton shower in the QGP in closed 
analytic form. Our result includes for the first time the quantum 
interference effects between the color currents associated with 
the primary parton and the radiated gluon. We demonstrated 
that the interference term that spoils the naive classical super- 
position Il2ll of incoherent energy and momentum transfer rates 
persists for time scales relevant for heavy ion phenomenology. 
We showed that this term vanishes in the t — > 00 limit, as it 
must. 

We presented two applications of the source term derived 
in this paper: the evaluation of the differential rate of energy 
transfer by a parton shower to the medium and the hydrody- 
namic response of the medium to a parton shower. The for- 
mer should not be confused with the rate of collisional energy 
loss IU4I1 of the parent parton itself, which remains constant 
and relatively small. We found that the rate of energy and and 
momentum transfer from the parton shower to the medium de- 
pends sensitively on the angle of emission and number of glu- 
ons emitted. For collinear bremsstrahlung, the medium can- 
not resolve the primary parton and the radiated gluon and there 
is pronounced coherent suppression in the magnitude of the 
multi-parton source term. For realistic large angles of emis- 
sion [17], we actually find an enhancement in its amplitude. 
Using the GLV formalism to describe the medium-induced par- 
ton shower 01811 . including bremsstrahlung gluon multiplicities, 
energies and emission angles, we estimated that for LHC appli- 
cations lEoll a 100 GeV parton can transmit more than 1/3 of 
its energy to the QGP. In the future, we plan to conduct more 
detailed numerical simulations that include realistic geometry 
and expansion of the medium and improved calculations of par- 
ton shower formation in the QGP 12511 . We note, however, that 
our findings suggest that the rate of jet energy dissipation in 
the medium may significantly complicate the experimental/jet 
medium background separation at RHIC and at the LHC ifl l^ql . 

We embeded the newly-derived source term in linearized 
hydrodynamics to study the response of the QGP to the in- 
medium parton shower. We found that for small gluon emis- 
sion angles and minimal shear viscosity to entropy density ratio 
rj/s = 0.08 12211 the parton shower generates a well defined 
Mach cone similar to what one would expect from a single pri- 
mary parton or from a phenomenological incoherent collinear 
superposition of a parton with gluons fll2ll. For realistic mul- 
tiple gluon emission at large angles IU7I1 there is an enhance- 
ment in the rate of energy transfer to the medium not only from 
the development of the parton shower but also from the inter- 
ference effects discussed in this Letter. This, in turn, yields 
a larger perturbation of the energy density in the medium but 
there is no longer a well defined Mach cone in the QGP. Our 
result is intuitive and easy to understand: each parton in the de- 
veloping shower gradually becomes a separate source of energy 
deposition and excitation of the strongly-interacting medium. 



The overlapping perturbations in very different spatial direc- 
tions wipe out any distinct Mach cone structure. Our finding 
is critical for heavy ion experiments and phenomenology. It 
suggests that schematic treatments of the parton shower pro- 
vide an inaccurate description of Shockwave phenomena in the 
QGP. It also suggests that high pj triggers are not likely to have 
well defined away-side associated Shockwave structures, since 
the number of medium-induced gluons grows with the parent 
parton energy 111 811 . 
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